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ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ 


УДК 539.375 
Б.И. СМЕТАНИН, Б.В. СОБОЛЬ, С.С. ВОЛКОВ 


ОБ ОДНОМ ЭФФЕКТИВНОМ МЕТОДЕ РЕШЕНИЯ 
СИНГУЛЯРНЫХ ИНТЕГРАЛЬНЫХ УРАВНЕНИИ ЗАДАЧ 
МЕХАНИКИ СО СМЕШАННЫМИ ГРАНИЧНЫМИ УСЛОВИЯМИ 


Рассматривается двумерное сингулярное интегральное уравнение в достаточно 
произвольной симметричной области, к которому сводятся некоторые задачи меха- 
ники со смешанными граничными условиями. К таковым, в частности, относятся 
задачи о плоском ударе пластинки о поверхность несжимаемой жидкости, а также 
задача о плоской трещине нормального разрыва в неограниченной упругой среде. 
Специальным представлением ядра интегрального уравнения удается существенно 
ослабить его сингулярность, что, в свою очередь, позволяет реализовать устойчи- 
вый вычислительный процесс по методу коллокаций. 

Ключевые слова: сингулярное интегральное уравнение, трещина, метод колло- 
каций, полином Чебышева П рода. 


Постановка задачи. Рассмотрим интегральное уравнение следующего 








вида: 
(<, п)В (у, 5,1) 

Ау =2лху; (ку) ЕО; в= Ев. @ 
лы Что бъз) ваз в е-уча-Ут. © 


Здесь © - некоторая замкнутая область, имеющая две оси симметрии Ох и 
Оу; Ч (& , 7) - неизвестная функция, удовлетворяющая условию 


4(х, у), =0. (2) 

В задаче о плоской трещине нормального разрыва уравнение (1) 

вытекает как следствие непосредственного интегрирования известного ин- 
тегродифференциального уравнения [1] с учетом симметрии задачи. 


В этом случае 9(Х, У) - приведенная к безразмерной нагрузке ^ 
амплитуда раскрытия трещины. Здесь р=соп<Ё- интенсивность нормальной 
нагрузки, приложенной к берегам трещины; 9 = Е[2(1-у’)] ', Е - мо- 
дуль Юнга, у - коэффициент Пуассона. 

Интегральное уравнение [1], в частности, также можно трактовать 
как модель задачи об ударе абсолютно жесткой плоской пластины формы в 
плане © о поверхность несжимаемой жидкости. В этом случае определе- 
нию подлежит функция импульсивного давления под пластиной. 

Если кривизна контура / области © , рассматриваемая как функция 


дуги $5, принадлежит Н!' (Г), а > 0, то решение уравнения (1) имеет вид 


а(х, у) = /Кх, у)®(х, у), (3) 


где [(х, у) = 0 - уравнение границы области ©. 
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В частности, в задаче о трещине размерная функция раскрытия 


ь ра 
связана с ней у(х, у) =^9 96, У), где а — постоянная, имеющая раз- 


мерность длины. 
В работах [2,3] предложена схема итерационного процесса по- 

строения последовательных приближений функции @(х, у). Слабым зве- 

ном этой схемы пока остается неустойчивая сходимость вычислительного 

процесса на осях симметрии. 

Метод решения. Учтем искомое решение (3) в (1) и представим его 

в виде 


[Е т п)К ит ‚1)(@(&,1)- @(х, у)) В 
—х)(п- у) 


о следующим представлением интегралов вида (1) в 
виде суммы двух Е 


ом (Ех)? + (п- у)? т=| (6, м) (6-х) аа 
= хи - ›) о(п- 8-х) +) 


(6,1) (т- У) Чат. (5) 
а(Е- ХЕ)? +(и- 5 
При таком представлении удается достичь под знаком интеграла 
снижение вдвое порядка особенности, а именно, разделение особенностей 
по &и по п; у каждого слагаемого остается только одна своя особенность. 
Приведем здесь также представление уравнения (1), рассматри- 
вавшегося в [4]. 


АЕ = лх». (4) 























@(х, Ур. +1 =2лху. (6) 
Здесь 
1 =0(х, ЭГ | к о ат ; 
рае ‚П я п) ®(х, у)) К ЧЕ 
ее -— а 
ы —х)(П — У) 
//(5,П ны 1) @(х,п))Е 
ат 
ы —х)(п- у) 


Выделение разностных множителей в данном случае будет в даль- 
нейшем использовано в технических целях при вычислении интегралов на 
одной из осей (для улучшения сходимости). 
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Рассмотрим интегральное уравнение (1) В области 
О:1-х? — у. > 0. Уравнение границы области соответственно имеет вид 


2 4 о 
Цх, у) =1-х’-у =0. Такая граница является внешней по отноше- 


нию к границе окружности /(х, у) =1-х* -—у*=0, совпадая с ней в 


вершинах (0,=1); (+1,0). Это обстоятельство в дальнейшем учитывается 
при реализации вычислительного процесса. 
В качестве нулевого приближения решения можем принять [2] 


@,(х, у), = 2лху. (7) 

При сделанных предположениях относительно области © функция 

4(х, у) является четной по обеим переменным. С учетом этого преобразу- 
ем уравнение (1) к виду 


[челн (&и»х, у)ават = 2лху. (8) 
Здесь о 


Н(Ельх,)= (о. РОК ОО, 


В,, = =)? +1); О, =(ЕХ у)”. (9) 
Нетрудно показать, что ядро интегрального уравнения (8) является 
четной функцией по & и 1 и нечетной по хи у. В результате преобразо- 
ваний нетрудно также убедиться, что в центре области и на осях симмет- 
рии выражения т в нулевом приближении решения соответственно имеют 


вид [2]: 








Же 0: ев | уИ(&,п) Ааёат ; (10) 


ху 2 (Е? +1?) . 


= 0<у<1: а" ИЕ 7%  атки) 
и- У += ау +=? 


узфочеч ОУ ‚п | 6х 6х феи (12) 

уе -х’ чт’ Ех +т 
а к непосредственному вычислению выражений (10)-(12). 
Функция, стоящая под интегралом в (11), имеет особенность при & =0. 


Для того чтобы избавиться от неё, проведём регуляризацию, а именно, 
прибавим и вычтем выражение, стоящее в этих скобках, при & =0.В 


результате имеем 
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В ты ПУ 
22 & |\/м-у’+8 
П+У 


а Ват + 
ПТУ 
Е" 


заитьвая В о что 


п [ет ав пи |] = ГЕ о. 


-а(т) = -с(5) 


где для © 1-х?-у>04@)=\1-П, а с()=41-°, и значения 


следующих интегралов: 


(13) 





[521(п — у) — %ет(п + УЛАЕам. 


[Дату - у) зи + УМ =47 


и и и- = —Л (по Адамару), (14) 


перепишем это а в виде 
| ИЕ | У 
ту +Е 


= — [517 (1 — у) — я вп(п + т + 2му. (15) 
И 


Реализация непосредственных вычислений в (13) показывает ус- 
тойчивый счет при достаточно больших значениях у. Далее для малых у 
(0О<у<0.1; 0.1 < у < [) преобразуем (13) к следующему виду: 


а р 9/1) 





зо ++ 
1 


(И- 75 + (1+ У)* ++) 57 + (п - У)" 
Исследование (16) с выделением точного значения для круговой 
области показало устойчивый счёт при 0<у< 0.1. 


Аналогично можно получить соотношения для оси у=0,0 <х<1, 
вытекающие из (12). 








а&ат. — (6) 
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А Ем) 
оо т? + (Ех)? м? + (&+х)? 
1 


, ЧЕ. 
(Е. + (&+х)? + (Е+х) ут? + (Е -х)? 


Численная реализация соотношения (17) показала устойчивый счёт 
для у=О0 и Ухе ©» . Перейдем к рассмотрению общего случая. 








$ |+ 


(17) 








В Г, помножим числитель и знаменатель на (Е -х)+(п-У)? 
и, разбив его на два интеграла, получим следующее представление: 


[ ‚= 
в © себе 

а(П- у) (8-х) +(-у) 
+ Ут 1-9" + МУ" де (18) 

а (Е-)\Е-х’ +’ 
Проведя в (18) регуляризацию при & = х и = У в первом и 
втором слагаемом соответственно и воспользовавшись равенством (14), 

получим формулу для х = 0 и у + 0 следующего вида: 


Щ_ № 2-м? [ с) +12) 
5 @-У |/&-5+м-у? 


зи — 2) ФУ" ЧЕ пу [ вт Фут + 


Е 4) 2 2 
ф 4 ()-5 (-У семи +2. (19) 
аа (&—х) Е->: +(п-Уу 


Таким образом, полученные соотношения в совокупности позволя- 
ют эффективно вычислять значения Г, для Ух, уЕО. Это, в свою оче- 


























редь, дает возможность определить @,(х‚ у). Значения @%(х,у) для 


центра области и осей приведены ниже. 
Для реализации метода коллокации представим искомую функцию 


в виде разложения по полиномам Чебышева П рода С,» (х): 


ву) => У Х„О„сди,с— те =. (20) 


Реализуем принятое представление в (7). В результате получим 
систему линейных алгебраических уравнений для определения коэффици- 


ентов разложения Х 


тп" 


пря бы >= 21х.у,;(ь)=0.№).  @0 


т=0 и=0 


593 


Физико-математические науки 








У; 
с(х,) 
Произведем далее соответствующие вычисления [> Г; в (6). Для 
этого в выражениях (7), (15), (16), (19) выделим разностные множители 
(@©(х,1)-— @(х, у)) и(®(Е,п)-@(х,1)) и воспользуемся представле- 
нием (20). В результате получаем систему линейных алгебраических урав- 





Здесь точки Х,, Е 0,1] - узлы Чебышева. 


нений для определения коэффициентов Хх р 
М 


М 
УХ, +Ь+Ь)=2у, — (@1=0.М. (22) 


т=0и=0 
Результаты вычислений для области <>: 1-х —у >0 представлены 
в табл.1-3. Расчеты проведены соответственно по 16, 25 и 36 узловым 
точкам. Столбцы @(х,у),, 4(х,у), соответствуют решению, полученному 
методом Лейдло [1]. В таблицах использованы — обозначения: 
[= = у=(-27)° (заметим соответствие значений парамет- 
а” 
ров х=1<>Е=1). 
Таблица 1 
Осевые значения искомых функций на осях х=0 





@(х, У). (ху) а» (х,у) Чо, у) @, (х, у) Ч» (х,у) 
и ТТ ТМТ1 





0 0 16 | 0.67039 0.67039 0.74934 0.74934 0.66904 0.66904 
25 | 0.67164 0.67164 0.74934 0.74934 0.66818 0.66818 
36 | 0.67138 0.67138 0.74951 0.74951 0.66945 0.66945 





0 0.25 | 16 | 0.65397 0.65269 0.73163 0.73020 0.65043 0.64916 
25 | 0.65571 0.65443 0.73182 0.73039 0.64963 0.64836 
36 | 0.65553 0.65424 0.73190 0.73047 0.65091 0.64964 





0 0.50 | 16 | 0.61427 0.59476 0.68560 0.66383 0.60351 0.58435 
25 | 0.61538 0.59584 0.68636 0.66456 0.60311 0.58396 
36 | 0.61568 0.59613 0.68635 0.66456 0.60406 0.58488 





0 0.75 | 16 | 0.56283 0.46535 0.62767 0.51896 0.54710 0.45234 
25 | 0.56450 0.46673 0.62787 0.51912 0.54533 0.45088 
36 | 0.56401 0.46632 0.62795 0.51919 0.54627 0.45165 





0 1 16 | 0.46180 0 0.56909 0 0.48597 0 
25 | 0.49863 0 0.57009 0 0.49562 0 
36 | 0.46141 0 0.57046 0 0.49640 0 





0.25 |0 16 | 0.67493 0.65350 0.75011 0.72629 0.67113 0.64982 
25 | 0.67318 0.65181 0.75004 0.72624 0.67007 0.64879 
36 | 0.67354 0.65215 0.74997 0.72615 0.67154 0.65021 





0.5 |0 16 | 0.69040 0.59790 0.75266 0.65180 0.67827 0.58740 
25 | 0.69079 0.59817 0.75254 0.65172 0.67708 0.58636 
36 | 0.69048 0.59797 0.75257 0.65174 0.67927 0.58826 





0.75 |0 16 | 0.72514 0.47963 0.75756 0.50108 0.69223 0.45786 
25 | 0.72679 0.48073 0.75757 0.50108 0.69024 0.45655 
36 | 0.72660 0.48060 0.75752 0.50105 0.69347 0.45869 








1 0 16 | 0.80227 0 0.76530 0 0.71312 0 
25 | 0.81346 0 0.76526 0 0.71716 0 
36 | 0.81400 0 0.76512 0 0.71875 0 
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Таблица 2 
Осевые значения искомых функций на осях у=0 
Хх |1 № | @(х,У), | аку), | 6% (У) | 9, (,У) | ®,(х,у)| Чк (ху) 
ТМТ1 М1 
0.25 | 0.25 | 16 | 0.65908 0.63691 0.73307 0.70841 0.65295 0.63098 
25 | 0.65765 0.63552 0.73310 0.70843 0.65198 0.63004 
36 | 0.65805 0.63591 0.73333 0.70865 0.65331 0.63133 
0.25 | 0.5 16 | 0.61775 0.57914 0.68831 0.64529 0.60622 0.56833 
25 | 0.61685 0.57829 0.68854 0.64550 0.60538 0.56754 
36 | 0.61724 0.57866 0.68794 0.64494 0.60657 0.56866 
0.25 | 0.75 | 16 | 0.56439 0.45182 0.63092 0.50507 0.54839 0.43901 
25 | 0.56415 0.45162 0.63046 0.50471 0.54731 0.43815 
36 | 0.56449 0.45189 0.63008 0.50441 0.54879 0.43933 
0.25 |1 16 | 0.50422 0 0.57271 0 0.48829 0 
25 | 0.51824 0 0.57190 0 0.48841 0 
36 | 0.49520 0 0.57149 0 0.48082 0 
0.5 0.25 | 16 | 0.67407 0.58262 0.73731 0.63728 0.66086 0.57120 
25 | 0.67310 0.58178 0.73734 0.63731 0.65932 0.56987 
36 | 0.67345 0.58209 0.73616 0.63629 0.66158 0.57182 
0.5 0.5 16 | 0.62841 0.52694 0.69603 0.58364 0.61552 0.51613 
25 | 0.62772 0.52636 0.69636 0.58391 0.61426 0.51507 
36 | 0.62764 0.52629 0.69570 0.58336 0.61863 0.51873 
0.5 0.75 | 16 | 0.57034 0.40837 0.64077 0.45881 0.55744 0.39914 
25 | 0.57044 0.40845 0.64152 0.45934 0.55856 0.39994 
36 | 0.57028 0.40834 0.64144 0.45929 0.55784 0.39942 
0.5 1 16 | 0.55096 0 0.58372 0 0.49352 0 
25 | 0.50933 0 0.58439 0 0.48531 0 
36 | 0.56303 0 0.58395 0 0.49224 0 
Таблица 3 
Осевые значения искомых функций вне координатных осей 
* { м @(х, У). а(сху) Г: @,(х,у) Чо(х, У) ©, (х, у) Чк (х, У) 
ТМТ1 ТМТ1 
0.75 | 0.25 | 16 | 0.70155 0.46312 0.74533 0.49202 0.67654 0.44661 
25 | 0.70337 0.46432 0.74546 0.49211 0.67392 0.44488 
36 | 0.70298 0.46407 0.74553 0.49216 0.67887 0.44815 
0.75 | 0.5 16 | 0.64565 0.41350 0.71175 0.45583 0.63671 0.40777 
25 | 0.64623 0.41386 0.71182 0.45587 0.63108 0.40416 
36 | 0.64552 0.41341 0.71213 0.45607 0.63513 0.40676 
0.75 | 0.75 | 16 | 0.58212 0.31835 0.66487 0.36360 0.58651 0.32074 
25 | 0.58407 0.31941 0.66473 0.36352 0.57428 0.31406 
36 | 0.58470 0.31976 0.66551 0.36395 0.58966 0.32247 
0.75 [1 16 | 0.48561 0 0.61375 0 0.52252 0 
25 | 0.47340 0 0.61377 0 0.61689 0 
36 | 0.47299 0 0.61496 0 0.47355 0 
1 1 16 | 0.53108 0 0.71964 0 0.69736 0 
25 | 0.55493 0 0.72921 0 -2.6032 0 
36 | 0.64114 0 0.73624 0 5.336 0 
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Аналогичным образом все полученные результаты можно распро- 


странить на область более общего вида о :1- (х/а)’ -(»/Б)* >0. При 
этом в процессе вычислений выделяется соответствующая эллиптическая 
область ©. :1- (х/а)? — (у/Б)?>0. 
Оценку интегральной характеристики проводили по формуле 
Ре [с у)@(х, у)ах4у (табл.4). В работе [5] приведена оценка ин- 
о 


тегральной характеристики для нашей области ( ©<>: 1 _У >0) Р<Р,, 
где Р, =1.565. 














Таблица 4 
Оценка интегральной характеристики 
ы Р, р 
16 1.54743 1.565 
25 1.54324 1.565 
36 1.55037 1.565 

















На рисунке представлен график искомой функции а(х,у) для об- 
ласти О: 1-х*-у>0. 





График функции а(х,у) 


Анализ численных результатов. 

1. Данные табл.1-3 свидетельствуют об устойчивой сходимости 
вычислительного процесса как на координатных осях, так и в других точ- 
ках рассматриваемой области. 

2. Данные табл.4 подтверждают достоверность полученных ре- 
зультатов ввиду выполнения необходимой оценки для известной инте- 
гральной характеристики. 

3. Полученное решение для задачи об ударе является весьма при- 
емлемым, но для задачи о трещине оно требует более высокой точности, 
так как в этом случае ключевым моментом является определение искомой 
функции на контуре с целью последующего определения коэффициента 
интенсивности нормальных напряжений. 
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